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Abstract 

We show that the supersymmetric extension of the quantum spherical model can be seen as a by-product 
of the stochastic quantization of the classical spherical model. This is so because of the well know equiva- 
lence between the Brownian motion described by means of the Langevin equation and the supersymmetric 
quantum mechanics, and it is connected with the existence of the Nicolai map. 
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I. INTRODUCTION 

The stochastic quantization procedure was proposed by Parisi and Wu as an alternative 

method to quantize field theories. Roughly speaking, the construction is based in an analogy 

_ s 

between Euclidean quantum field theory whose functional integration is weighted by the factor e R , 

with S the Euclidean action, and the Boltzmann distribution of a statistical system in equilibrium 

u 

e , where U is the energy. By the proposal, to the quantum field it is associated a stochastic 
process whose evolution follows a Langevin or Fokker-Planck equation. The Green functions of 
the quantum theory are then obtained as the equilibrium limit of the correlations functions of the 
stochastic theory. Various investigations have been reported on its conceptual/technical aspects 
and the method has been applied to numerical simulations in lattice field theories (see 0-13 and 
the references therein). 

From another point of view, stochastic quantization provides connections between field theoret- 
ical models. Indeed, it relates a classical theory in D dimensions with a theory in Z) + 1 dimensions, 
because it introduces an additional "time" direction in the system, a fictitious time whose evolution 
is governed by the Langevin equation. In this sense, stochastic quantization would be interpreted 
as a prescription for the construction of higher dimensional field theories, without taking the equi- 
librium limit and consequently considering the fictitious time as a physical parameter. A modern 
and related approach of these aspects is presented in Qj. 

The last observations may be positively exploited in favor of the quantization of the classical 
spherical model Q], as explained from now on. The Hamiltonian of the model, in the presence of 
an external field h, is given by 



-'^^Jr,r'SrSr'+h'^^Sr, (1) 



where r and r' are lattice vectors, {S^} is a set of spin variables that can assume continuous values, 
—oo < Sr < oo, in a D-dimensional hypercubic lattice with periodic boundary conditions; 
is the interaction energy that depends only on the distance between the sites r and r', Jj.y = 
J(|r — r'l). The variables are subject to the spherical constraint, 

T.Sr=N, (2) 
r 

where N is the total number of lattice sites. The first step towards the quantization of the model is 
to provide a dynamic to the system. This can be done by including in the above Hamiltonian a term 
involving the conjugated momenta of the spin variable or, equivalently, by adding to its Lagrangian 



a kinetic term containing time derivatives of the spin variables. However, there is no unique way to 
choose the form of these terms and, moreover, different forms yield different dynamical behavior ^- 
Q]. The stochastic quantization furnishes a natural framework to ameliorate this situation, because 
the spin variables are assumed to satisfy the Langevin equation, which then gives rise to a dynamic 
for the system, with the identification of the fictitious time as the ordinary time. As a result of this 
process we are led naturally to a super symmetric theory. Complementing our recent work on the 
supersymmetric spherical model [icj and proceeding with applications of field-theoretical methods 
to the study of statistical mechanics, we intend here to discuss in detail these issues. 

Our work is organized as follows. We start by reviewing some aspects of stochastic quantization 
in quantum field theory and in quantum mechanics, that are important for discussions in the sub- 
sequent sections. In section III, we study the stochastic quantization of the mean spherical model 
and compare it with other approaches. Section IV is devoted to the construction of supersymmetric 
quantum spherical model from the classical one. A summary of the work is presented in section V. 



II. SOME REMARKS ON STOCHASTIC QUANTIZATION 

Before considering the stochastic quantization studies of the spherical model, let us initially 
outline some basic steps of the prescription for the stochastic quantization of a scalar field. Then 
we discuss the equivalence between Langevin equation and supersymmetric quantum mechanics. 
This simple analysis is an important point of our discussions to be performed in the next sections. 



A. Stochastic Quantization of a Scalar Field Theory 

In a general setting, the stochastic quantization of a field theory model, with a scalar field (pix) 
depending on D-dimensional Euclidean space-time coordinates x and specified by an action S, 
proceeds as follows. The field </)(x) is assumed to depend on a fictitious time r, (p{x) — ?> (/)(x,r), 
such that the evolution in this variable is governed by the Langevin equation 

(x,t) k 6S 



where n is the diff^usion coefficient that can be conveniently chosen (to facilitate the comparison 
with other works), r/(x,r) is a noise with a Gaussian distribution, 

(r/(x,r))=0, {ri{x,T)ri{x',T')) = kS{t-t')6^{x-x') (4) 
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and S is the Euclidean action that must include also an integration over the fictitious time, 

S = J d^xdT£{(P{x,T),df,(f){x,T)). (5) 

The theory arising from this process is a D + 1-dimensional one. Let us describe how to come back 
to a D-dimensional but quantum theory. For a specified initial condition in the fictitious time, the 
Langevin equation ([3]) has the noise-dependent solution 0^(x,t) whose correlation functions are 
defined as 



t>n{xi,Ti) ■ ■ ■ (j)r,ixN,TN))ri = j T^V (pr/ixi, Ti) ■ ■ ■ (j)r,{xN , Tn)P[ii], 



with the normalized probability 



= — exp 



/ d^xdr rf'{x,T) 

2k J 



(6) 



(7) 



and partition function 



Z = I Vrj exp 



1 

'2k 



d^xdr ri^{x, r) 



(8) 



The equivalence of this procedure with other quantization methods is established at the level of 



the Green functions 



In fact, it has been shown that under certain conditions which include 
the semi-positivity of the associated Fokker-Planck Hamiltonian, if we take ri = • • • = rjv = r 
and then take the equilibrium limit r — > oo, we obtain the time-ordered Green functions of the 
quantized field theory, 

lim ((/.(xi, r) • • • (/>(x„, r))^ = (0|T(/.(xi) • • • 0(xjv)|O). (9) 

r— ^oo 

Shortly, the quantization procedure above described can be summarized as follows. We start 
with an action in D-dimensional Euclidean space-time. The field is assumed to depend on a new 
time coordinate so that the resulting theory is a D + 1-dimensional one. In order to reach the 
quantum theory we take the equilibrium limit so eliminating the fictitious time and returning to a 
D-dimensional situation. 

However, as argued before, we may take a different route and consider this process as a prescrip- 
tion for the construction of D -|- 1-dimensional theories from D-dimensional ones. In this work we 
are interested in the D + 1-dimensional emerging theory itself. In quantum mechanics, for example, 
this establishes a quite interesting and well know connection between the Brownian motion and 
the supersymmetric quantum mechanics, to be presented now. 
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B. Langevin Equation and Supersymmetric Quantum Mechanics 



Let us start by considering a point particle whose spatial coordinate is x. Following the stochas- 
tic quantization prescription, the coordinate x is assumed to depend on the fictitious time r, 
X — > x{t). As mentioned before, we will not eliminate the fictitious time in the end of the process, 
but instead we will identify it just as the ordinary time. The Langevin equation (with k = 1) reads 



dx{T) 
dT 



16S , , 



(10) 



that describes a Brownian motion of the particle in a heat bath subjected to a random noise. The 
partition function is given by the integral over the noise 



Z = Dr/exp 



-\ I dr v'{r) 



(11) 



We perform a change of variable of integration r/(r) — )• x(r) defined by the Langevin equation p!Ol 
to write the partition function as 



"Dx det 



6x{t') 



(12) 



The determinant can be written as an integral over anticommuting (Grassmann) variables '0 and 



det 



6x{t') 



'DipT)'il) exp 



(13) 



where we have introduced ^ = ^ff with V = With this, the partition function becomes 



'DxVil^'Dip exp 



dr 



\ f dx\ .1^,2 T ( d 



(14) 



discarding surface terms. From this expression, we may read out the effective Euclidean Lagrangian 

(15) 



that is exactly the Lagrangian of the supersymmetric quantum mechanics |lll.[l^. As can be easily 
verified, it is invariant under the supersymmetry transformations 



5eX = —tpe, Sf:tp = x e — V e and d^tp = 



(16) 



and 



6^x = eip, 6eip = and 6^^jJ = e x + eV. 



(17) 
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As usual, this super symmetry can be cast in a off-shell formulation by introducing an auxiliary 
field F, so that the Lagrangian becomes 

From the off-shell supersymmetry it is straightforward to construct the superspace formalism, which 
for the quantum mechanics it is constituted by the anticommuting coordinates 6 and 9 besides the 
time T. For example, the superfield is expanded in terms of 9 and 9 

<^ = x + ip9 + 9ip + 99F. (19) 

With the supercovariant derivatives 

we may write the superspace action from the Lagrangian 

C = ^D<i>D<^> + L{<i>) =^ S = jdTded9C, (21) 

such that L' = V. In general, the existence of local or global symmetries implies relations between 
the Green functions of the corresponding theory, designated generically as Ward identities. In 
particular, that is the case for the supersymmetry. A natural question then is about the role of the 
supersymmetry in the case discussed here. Interesting enough is the fact that this supersymmetry 
is the reflex of the forward and backward time propagation symmetry in an equivalent Fokker- 
Planck formulation, corresponding to the occurrence of two Fokker-Planck Hamiltonians known in 
the literature as forward and backward Hamiltonians [l^. Nevertheless, in the equilibrium limit 
r — > oo only the Hamiltonian associated with forward time propagation will be relevant. 

The procedure outlined above is related with the existence of the so called Nicolai map [14[. The 



Nicolai map deals with the characterization of super symmetric theories by means of functional inte- 
gration measures. Roughly speaking, if we start with a super symmetric theory and then integrate 
out the fermionic fields from the path integral, this will produce a non-trivial determinant. There 
exists a transformation (the Nicolai map) for the bosonic fields whose Jacobian cancels exactly 
the determinant coming from the fermionic integration. The result is just a bosonic free action. 
In the quantum mechanical case, the Langevin equation constitutes an explicit realization of the 
Nicolai map! The route from Langevin equation towards the supersymmetric quantum mechanics 
corresponds to a Nicolai mapping in the inverse direction. 
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III. STOCHASTIC QUANTIZATION OF THE MEAN SPHERICAL MODEL 



In this section, we wish to discuss the stochastic quantization of the mean spherical model 
whose constraint is imposed thermal average 

J2r{Sr) = N (not strictly as in ([2])) in contrast 
with other approaches like the canonical and the path integral quantization methods applied to 
this specific model 3: S fisl 1 



17|. Let us consider the classical Hamiltonian 



r r,r' \ r / 



(22) 



where 5 is a positive parameter associated with the dynamics. In fact, the presence of the kinetic 
term is only to provide a dynamic to the model, but we will not impose commutation relations 
between and P^- The quantum character comes from the stochastic process. The corresponding 
Lagrangian in the Euclidean time (which includes a global signal redefinition) is given by 



= ^ E ^rit) + ^ E Jr,r'Sr{t)Sr>{t) + /i ( E (0 " n] 
^ r r,r' V r / 



(23) 



We shall use the variable t to represent the Euclidean time and the variable r will be reserved 
to the fictitious time in the stochastic quantization process. We start the stochastic quantization 
procedure by assuming that the spin variables depend on the fictitious time Sr{t) — )• Sr{t, r), whose 
evolution is governed by the Langevin equation (with k = 2), 

dSr{t,T) _ 5S 



dT 



SSr{t,T] 



(24) 



where S is the Euclidean action. 



S 



J ^E'^r(*'^) + ^E^r,r'Sr(t,r)Sr.(i,r) + M^E'^r(i,^)-A^^ 



(25) 



and rij.{t,T) is the Gaussian noise. It is worth to stress that the Langevin dynamics considered 
here constitutes just an intermediate step towards the quantization of the model. Notice that 
the Langeviii dynamics itself for the classical spherical model was studied by several authors, for 
example [18h21|. 



We will consider the problem of stochastic quantization at finite temperature, that can be 
put in consonance with the imaginary time formalism approach of field theory. In this situation, 
the Euclidean or imaginary time is restricted to the interval [0, /?] and the fields must satisfy 
periodic or anti-periodic conditions according to their bosonic or fermionic character, respectively, 
i.e, (p{0) = ±ip{/3). Besides the fields, the noise must also satisfy periodic or anti-periodic conditions. 



Namely, for a bosonic or fermionic noise we have r/(0) = ibr/(/3), assigning to it the same character 
as the corresponding field. Then, the Gaussian noise average must be adequate to reflect the 
periodic conditions of the bosonic fields in the imaginary time i.e.. 



(r?r(t,r))=0 and {,^,(t,T)rjr'{t' ,t')) = - ^"'"^''''^SrySir - r') 



(26) 



where a;„ is the Matsubara frequency cvn = 27rn//3, with n € Z. It is convenient to introduce 
the Fourier decomposition of the variables Sr{t,T) and r]j.{t,T) taking into account the periodic 
conditions. 



Sr{t,T) 



^g.(r.qW)5^(^^^^) and Vr{t,T) 



Y^S'-^+^-'\^{uJn:T). (27) 



q.n ^ ^ q,n 

By identifying J(q) as the Fourier transformation of the interaction energy J,. ^ 



J(q) = J;J(|h|)e-^•^ 



with h = r — r', the Langevin equation in the Fourier space takes the form 



(28) 



g5'q(c<Jn,r) 

dT 



+ 



+ 2^ + J(q) 



S'q(w„, r) - 7/q(a;„, r) = 0, 



whose general solution is given by 

S^{un, r) = S^icon, 0) e-^'^'i'")^ + / df e-f^'('i.")(---) r?q(u;„, f) 
where J7(q, n) is defined as 

J72(q,n) = 



(29) 



(30) 



(31) 



From the solution (I30p we can determine the correlation functions of the mean spherical model 
according to (|26p . As stated early, the correlation functions so obtained must recover those of 
the quantized spherical model whenever we take the appropriate limit of the n-point function 
Ti = r2 = • • • = Tn = r — > oo. In particular, for the 2-point function we get 



1 



-n2(q,n)|T 



(32) 



In this expression we have discarded terms that do not contribute in the equilibrium limit. Then, 
by taking t = t' oo, we find the correlation function in the momentum space for the quantum 
spherical model 

1 g 



(5q(w„)S'q/(w^)) = 5c^-c^'Sn, 



/32 [ul + 



(33) 



where 



2 

Consequently, the correlation function in the coordinate space is given by 



wl^2g{,.+ ^My (34) 



q ^ 

The equal-time correlation function reduces to, 

that is the usual quantum pair correlation function of the spherical model (see, for example jl6l|). 
We can also obtain the classical correlation function by taking the limit g( ^ 0: 



2 J 

The quantum self-correlation furnishes the constraint condition corresponding to the mean spherical 
model if we impose X]r('^r) = ^ i 

1-lV^cothf^Uo. (38) 

The constraint equation (|38p may be used to determine the critical behavior of the model 
To this purpose it is convenient to parametrize the interaction energy as J(q) ~ |q|^ for small 
values of |q|. The parameter x specifies the range of interactions, for example, for short-range 
interactions we have x = 2. The quantized spherical model is quite interesting in the sense that 
it exhibits a quantum phase transition, e. g., a phase transition at zero temperature besides the 
usual finite temperature phase transition. 

To sum up, in this section we applied the stochastic method as an alternative way to quantize 
the spherical model showing how to obtain the quantum pair correlation function as well as the 
constraint condition after the elimination of the fictitious time. Before closing, we sketch the 
procedures of canonical and path integral quantization of the mean spherical model discussed by 
several authors 



171 ]. In the canonical quantization scheme, the variables Sr and Pr of the 
Hamiltonian (|22p are promoted to operators analogous that of position and momentum, satisfying 
the commutation relations 

[5r, -Pr'] = i-^r.r' and [^r, ^r'] = [-Pr, -Pr'] = 0. (39) 
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The remaining task is the diagonaUzation of the harmonic oscillator-type Hamiltonian by means of 
the Bogoliubov transformations. The path integral formulation in the imaginary time formalism, 
on the other hand, involves a functional integration over the spin variables 

Z = f YlVSrC-^o^, (40) 

r 

with C the Euclidean Lagrangian (j23p . Using a Fourier decomposition for Sr{t) similarly to (j27p . 
taking into account the periodic condition S'r(O) = Sr{/3) required for bosonic variables at finite 
temperature, it is straightforward to perform the Gaussian integrals and determine the partition 
function. From either of these methods we may find the pair correlation function (|36p and the 
constraint condition ([38|l . 



IV. SUPERSYMMETRIC SPHERICAL MODEL 



Based on the existence of a Nicolai mapping, in this section we will obtain the supersymmetric 
quantum version of the spherical model from the stochastic quantization procedure of its classical 
counterpart. Different of previous case, here we will not consider a dynamical model from the 
beginning. In fact, as we mentioned in the Introduction, there is no necessity of the addition of a 
kinetic term. On the contrary, the dynamic arises from stochastic quantization scheme as long as 
we recognize the fictitious time as the ordinary time. For the strict spherical model, to be discussed 
here, we need to implement the stochastic quantization for constrained systems 



22 



23[, in which 



the Langevin equation itself is considered as a constraint. Our approach is similar to that for the 
case of the field theoretical non-linear o"- model {24 ] . 

Let us start with the classical Hamiltonian of the spherical model without external field 

'Hc = Y,KrySrSr', (41) 
r,r' 

with an interaction energy designated here by Kj-y = i^(|r — r'|) and under the spherical constraint 

Y,S'r=N. (42) 

r 

The spin variables are assumed to depend on the fictitious time that we designate by t, Sr — > Sr{t), 
and the action is given by 



S 



J dtC = j dt Kr,r'SrSr' + a ~ ' '"^^^ 
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In the functional approach, by integrating over the auxihary field a means that we are dealing with 
the strict spherical model and not the mean spherical model. The Langevin equation is written as 

dS^ 1 5S 



(44) 



dt 2 

where rjr is the Gaussian noise of each lattice site and we are using now n = 1. After calculating 
the functional derivative of the action with respect to the Sr variable, the Langevin equation takes 
the form 



- ^ Kr,r'Sr' - aSr + ?7r- 



(45) 



It is convenient to introduce a specific notation to unify the Langevin and the constraint equations. 
To this end, first we consider the set of N + 1 fields (S'r, cr) as components of a variable Sa- 



Sr, a = l,2,--- ,N 
a, a = N + 1 

Next, we introduce the functions 



r)r, a = l,2,--- ,N 
0, a = N + l 



(46) 



FriSr,a), a = l,2.---,N 
F{Sr), a = N + l 



(47) 



with the definitions 



and 



Fr{Sr, a) = Sr+J2 ^r,r"5r' + aSr 



F(5r)^^52-iV. 



(48) 



(49) 



According to this notation, the Langevin and the constraint equations are unified into the equation 
J^((S) = M. The partition function is given by the integral over the noise 



(50) 



The integration measure VJ\f symbolically stands for the functional integration over all the site 
lattices, VAf = Yl^^-^a — Tlr^^^- We will adopt a procedure similar to that of Faddeev-Popov 
path integral quantization for gauge theories. With the goal of making a change of integration 
variable from J\f to S according to ^(«S) = M, the partition function can be written as 



(51) 
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Prom a fermionic integral representation of the Faddeev-Popov-like determinant, 



det 



J {V^){V^) exp J dtdt' J2 ^ 

a.b 



(52) 



with 



and '^a = < 



^C, a = N + l 
and the integral representation of the delta function, 



V^r, a= 1,2,--- ,N 
C, a = N + l 



(53) 



(54) 



with 



Ar, a = 1,2, ••• ,Ar 
a, a = Ar + l, 



(55) 



the partition function acquires the form 

Z = j {VJ\f){VS){Vk){V^){V^ 



X exp < 



+ y dtdt'^^.^v&A. (56) 
a,6 I 



By integrating over the variable M (Gaussian) and using the fact that Ma = for a = iV + 1, we 

have 



Z = I {VS){VK){V'^){V^) exp 



(57) 



a=l a a,b ^ 

This result is the master functional and when written in terms of the components it becomes 

Z = J {VS){VX){Va){Va){VtP){Vij){VC)(VC) e-^'^'^'*, (58) 



with the master action 



^mast 



r r r,r' r \ r / 



(59) 
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Next, by integrating over the fields a, a, C, and C in equation (|58j) we end up with 



X exp 



• (60) 



In order to cast this expression in the desired final form, we need to perform a shift in the auxiliary 
field in the functional integration Ar — > Ar + (that does not affect the constraint because 
SrSr = f ^(X^r "^r) = 0) and also it is convenient to change t — )• — t in the action. After this, 
we get the partition function 

y r r r,r' r r,r' J 



X exp 



■ (61) 



This expression is interesting due to the following reasons. First, we may read out the Euclidean 
Lagrangian 



(62) 



r,r' 



that is invariant under the set of supersymmetry transformations 



and 



5eSr = V'rC 'JeV'r = SrE + Are, (J^Vr = and ^e^r = " V'rC 



5iSr = eipr, Silpr = 0, S^lpr = —SrS + Xr^ and 5eXr = ^Ipr 



(63) 



(64) 



where e and e are fermionic infinitesimal parameters. Moreover, the delta functions in the integrand 
impose the constraints 



(65) 



S'r = A^, Y^ = 0, V'r'S'r = 0, and Ar^r = V'rV'r 



that turn out to be the supersymmetric constraint structure which are closed under (|63p and 
From these observations, we identify the model above obtained as an Euclidean off-shell version 
of the supersymmetric extension of the spherical model presented in {lo| . In the later case, a 
supersymmetric version was constructed directly by adding fermionic degrees of freedom at each 
lattice site besides the bosonic ones with the purpose of studying its critical behavior at finite as 
well as zero temperature. 
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Here, in its turn, the super symmetric version of the quantum spherical model arises from the 
classical model by means of the stochastic quantization prescription. This comes from the con- 
struction of a D + 1-dimensional theory from a D-dimensional one adapted to constrained systems 
and with the properly identification of the fictitious time as the ordinary time. Interesting enough 
it is the arising of the precise super symmetric structure of the constraints. Finally, it is worth 
to say that this structure can straightforwardly be put in a superspace formulation, as presented 



m 



V. SUMMARY 

In this work, we reported some results about the stochastic quantization of the spherical model. 
We started by reviewing some basic aspects of this method with emphasis in the connection between 
the Langevin equation and the supersymmetric quantum mechanics, aiming at the application of 
the corresponding connection to the spherical model. An intuitive idea is that when applied to 
the spherical model this gives rise to a supersymmetric version that is identified with one studied 



m [10|. Before investigating in detail this aspect, we studied the stochastic quantization of the 
mean spherical model that is simpler to implement than the one with the strict constraint. We 
also highlight some points concerning more traditional methods discussed in the literature like 
canonical and path integral quantization. 

To produce a supersymmetric version, grounded in the Nicolai map, we investigated the stochas- 
tic quantization of the strict spherical model. We showed in fact that the result of this process is an 
off-shell supersymmetric extension of the quantum spherical model (with the precise supersymmet- 
ric constraint structure). That analysis establishes a connection between the classical model and 
its supersymmetric quantum counterpart. The supersymmetric version in this way constructed is 
a more natural one and gives further support and motivations to investigate similar connections in 
other models of the literature. 
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